Three-dimensional quasi- Tonks gas in a harmonic trap 
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We analyze the macroscopic dynamics of a Bose gas in a harmonic trap with a superimposed two- 
dimensional optical lattice, assuming a weak coupling between different lattice sites. We consider 
the situation in which the local chemical potential at each lattice site can be considered as that 
provided by the Lieb-Liniger solution. Due to the weak coupling between sites and the form of the 
chemical potential, the three-dimensional ground-state density profile and the excitation spectrum 
acquire remarkable properties different from both fD and 3D gases. We call this system a quasi- 
Tonks gas. We discuss the range of applicability of this regime, as well as realistic experimental 
situations where it can be observed. 



Recent developments in low-dimensional trapping 
loading of optical lattices [^|^, and external modification 
of the interparticle interactions j3| , open fascinating per- 
spectives towards the achievement of strongly-interacting 
systems in cold atomic gases. In this sense, a remarkable 
example is provided by the recent observation of the su- 
perfluid to Mott-InsulatorfMI) transition in cold bosonic 
gases in optical lattices Other strongly-interacting 
systems in cold atomic gases have been recently consid- 
ered theoretically, as the case of large scattering length 
[3, rapidly rotating Bose gases H , and one-dimensional 
systems, both bosonic and fermionic p^ . 

One-dimensional systems differ significantly when 
compared to 3D ones. Most remarkably, in ID gases the 
interactions become relatively more important the more 
dilute the gas is. In particular, for very low densities or 
a large value of the s-wave scattering length a, the sys- 
tem enters the so-called Tonks-Girardeau (TG) regime, 
in which the bosons acquire fermionic properties 13, 14]. 
For finite interactions, the problem of delta-interacting 
ID homogeneous Bose gases was exactly solved by Lieb 
and Liniger (LL) by means of Bethe Ansatz . A com- 
bination of the LL solution, the local density approxi- 
mation, and hydrodynamic equations have been recently 
employed in the analysis of different properties of trapped 
ID Bose the ground-state density [7|, the exci- 

tation spectrum the non self-similar expansion in a 
ID guid e 0, and the correlation and coherence proper- 
ties [lil. The latter has also been recently analyzed by 
means of quantum Monte Carlo methods lllf. 

The accomplishment of strongly-interacting ID gases 
requires tight transversal trapping, low atom numbers, 
and possibly the increase of a via Feshbach resonances 
[j, 16, 17]. In this sense, 2D optical lattices are specially 
favorable, since the on-site transversal confinement can 
be very strong. For sufficiently intense lattices, the tun- 
neling rate could be made small enough to ensure that 
each lattice site behaves as an independent ID experi- 
ment. This paper analyzes the situation in which the 
tunneling is non negligible, but sufficiently weak to as- 
sume the ID character of the local chemical potential. 
In this regime, which we call the quasi- Tonks regime, we 



predict that the combination of tunneling and ID local 
chemical potential results in important modifications of 
the 3D ground-state density and the excitations. 

In the following we consider a gas of N bosons in a 
cylindrically-symmetric harmonic trap, with axial (ra- 
dial) frequency ujz We assume a superimposed 
periodic potential provided by an optical lattice of the 
form V\{x) 4- V\{y) — Vo{s\t? qx + sit? qy), where q — 
27r/A is related to the laser wavelength A, and Vb is 
the lattice amplitude. The lattice has a periodicity 
d = 7r/g — A/2. From a macroscopic point of view 
it is convenient to introduce the effective mass, m* , 
to describe the effects of the lattice on the dynamics. 
The value of m* can be obtained from d'^Eo/dk^ = 
h?/m*, where Eo{k) is the dispersion law corresponding 
to the lowest energy band. In tight-binding regime, m* 
can be related with the tunneling rate J, as m/m* — 
tt'^J/Eh, where Er — h^q^ /2m is the recoil energy, 
J = — / Wi{x){—h^Vl/2m+V\{x))wi+i{x)dx, and {w^} 
are the Wannier functions for the lowest band. 

We first consider the case of a single isolated lattice 
site, which can be well approximated as a dilute gas of 
Nt bosons confined in a very elongated harmonic trap 
with radial and axial frequencies uip and ojz {ojp ^ oJz)- 
The transversal confinement is strong enough to fulfill 
that the interaction energy per particle is smaller than 
the zero-point energy hujp of the transversal trap. In 
this way, the transversal dynamics is effectively "frozen" 
and the system can be considered ID. The Hamiltonian 
describing the physics of the delta-interacting ID gas is 
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is the Hamiltonian in absence of the harmonic trap, and 
giD — —2h^/raaiY). The scattering problem under one- 
dimensional constraints was analyzed in detail by 01- 
shanii 111, and it is characterized by the ID scattering 
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length fliD = (— a^/2a)[l ~ C{a/ ap)], with a the 3D scat- 
tering length, Op = yj2h/mujp the oscillator length in 
the radial direction, and C = 1.4603 .... For the ther- 
modynamic limit, a ID gas at zero temperature with lin- 
ear density n, is characterized by an energy per particle, 
which can be obtained from the corresponding LL inte- 
gral equations 0|. Assuming that the density variates 
sufficiently slowly, at each point z the gas can be consid- 
ered in local equilibrium, and the local energy per particle 
is provided by the LL equations for the density n{z). We 
call this approach Local Lieb-Liniger (LLL) model. 

Although, strictly speaking, the LLL approach is only 
valid for ID systems, for sufficiently low J, the virtual 
processes involving transitions into nearest sites, should 
just result in a small correction to the ID scattering of the 
order of J//i, where /i is the chemical potential. There- 
fore, if J//i ^ 1, the actual local chemical potential can 
be well approximated by that obtained from the LLL ap- 
proach. If the tunneling becomes very small, even for 
large values of Nt the system could enter into the MI 
regime. However, this regime typically demands an ex- 
tremely small value of J//i. Summarizing, in what we 
called the quasi- Tonks regime, l/iVj < J / fji <^ 1, the LLL 
chemical potential can be employed, and at the same time 
the tunneling cannot be neglected. 

Assuming a sufficiently slow variation of the gas den- 
sity, the dynamics can be well described by means of the 
corresponding macroscopic hydrodynamic equations |l8l | : 
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where n{x, y, z) is the macroscopic 3D density and v = 
[vx, Vy,Vz) is the velocity field. Note that all the informa- 
tion about the lattice is contained in m* , and in the local 
chemical potential n{n). The use of the hydrodynamic 
equations is j ustified by the superfluidity of the quasi- 
ID tubes [131 and by the avoidance of the MI phase. In 
the equation ^ the local chemical potential of the 3D 
system is related to the ID chemical potential, hid, pro- 
vided by the LLL approach. In order to obtain fiiD the 
value of aiD must be obtained for each site. The latter 
demands the knowledge of the effective oscillator length 



flp associated with the on-site radial confinement. Since 
the latter can be very large, a Gaussian Ansatz can be 
assumed for the on-site radial wave function. In this way, 
we can calculate Op by minimizing the radial energy 
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In Eq. IjlOl) we have neglected the interaction energy, 
which for Vq/Er > 3 does not provide any significant 
contribution. Once known the value of aio, we evalu- 
ate fj,iD{n) by numerically solving the corresponding LL 
integral equations. 
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Figure 1: Density profile as a function of the dimensionless 
radius — ^f{ri/Ri)^ (see text). The solid, doted and 
dashed lines correspond to log A — —3, 0, 3, respectively. The 
resulting fit with a power law dependence (1 — f'^)" provides 
s — 0.54, 0.72, 0.91, respectively, being an indicator of the 
transition from the MF regime to the strongly-interacting one. 

By imposing the equilibrium conditions dn/dt — Vi = 
0, we obtain the equation of state 



^j,[no{x, y, z)] ^ ^J,T - V{x, y, z), 



(11) 



where tiq is the equilibrium density, and fix is the chem- 
ical potential of the system. Inverting Ullfl . one obtains 
the expression for hq. Imposing the normalization to the 
number of particles N we obtain: 
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= 4<3 f ff,-^[^^^i-e)]dt, (12) 
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where ^ = 2/xm|aiDp//i^. From Eq. Ill 211 we observe 
that, similarly to the strict ID case 0,B^3i ground 
state of the system is completely characterized by a sin- 
gle parameter A. For ^ ^ 1 the mean-field (MF) regime 
is retrieved, whereas the TG is found for A <^ \. In 
Eq. (fT^ . jliD{x) = 277i|aiDp/^i£i/fi.^, is a function of 
the dimensionless density 77 ~ \aiu\d?"n- Fig. ^ shows 
the density as a function of the dimensionless radius f, 
where = / Ri)'^ , with Ri the Thomas- Fermi radii. 
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In order to compare the different cases, we re-scale the 
density with respect to the central one. Nevertheless, 
we stress that the values of Ri and the central density 
depend on the parameters, and in particular on the in- 
teraction regime. 
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Figure 2: Parameter / as a function of log A (see text). 

From the previous discussion, it becomes clear that the 
ID properties of the local chemical potential are trans- 
ferred to the 3D stationary shape of the cloud. In partic- 
ular, for small values of A, the transversal density profile 
is not an inverted parabola, as expected for the case of a 
BEC in a harmonic trap. We stress once more, that this 
remarkable effect is just possible in the presence of the 
2D lattice, which guarantees, under the above discussed 
conditions, a LLL chemical potential. 
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In the expression for the frequencies we employ the pa- 
rameter / = ^"(I), which is an universal function of the 
parameter A. We show this dependence in Fig. |21 The 
parameter / ranges from 2, at the TG limit, to I, for 
the MF regime. In this two limiting cases the function F 
can be obtained analytically, being F = p"^ (TG) and 
F = (MF). 

The quasi Tonks regime is not only achievable for real- 
istic conditions, but, actually, it is expected to be the case 
for typical parameters in ongoing experiments, as those 
of Fig. 3, which shows the dependence of the frequen- 
cies ()17I19|) on the lattice amplitude Vq/Eh- We stress 



In the following we calculate the frequencies for the 
lowest excitations of a quasi- Tonks gas. In particular, we 
evaluate the breathing and quadrupole modes, by using 
the Ansatz n = no(ri/bi) / Y\jbj, Vi = (bi/bi)ri, which 
constitutes an exact solution of the continuity equation 
(PJ. Multiplying the Euler equations by nxi and inte- 
grating, one obtains 

b,+Culb, + '^Fi Yl bj)^0 (13) 

by+^lby+'^Fi Yl b,)^0 (14) 

h+Lulb^ + '^Fi Yl b,) = 0, (15) 

where we have employed the effective frequencies u)\ ~ 
Lo\rn./rn* and defined the function 

^^^^ = /"oX.^Mp-^-o)dV. (16) 

Using the previously calculated equilibrium solution, one 
obtains that F{p) is independent of the particular choice 
of the coordinate i, and that F{1) = —1. Note that the 
function F only depends on A, which, as previously dis- 
cussed, completely characterizes the static properties of 
the system. Linearizing around the equilibrium solution 
bi = \, one obtains the frequencies of the breathing mode 
{liJb)i the m — Q quadrupole mode {i^qi) and the m = 2 
quadrupole one {u!q2)- 



(17) 

(18) 
(19) 

that the system enters the MI regime for large values of 
Vq/Ej^ (around 40 for the case of Fig. 3). If this is the 
case, as discussed above, our macroscopic hydrodynamic 
approach fails, and the system becomes a set of indepen- 
dent ID gases. On the other hand, for decreasing values 
of Vo/En the system abandons the quasi- Tonks regime 
(which for the case of Fig.Ooccurs at Vq/Eh ~ 15), and 
the condition J/ fi ~ 1 is reached. This constitutes an ad- 
ditional intermediate cross-dimensional regime, in which 
the gas is 3D, but the local chemical potential is not that 
expected for a 3D Bose gas. For this regime, contrary to 
the case of the quasi- Tonks gas, the local chemical poten- 
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Figure 3: Frequencies of the breathing mode (solid), 
quadrupole 1 (dotted) and quadrupole 2 (dashed) as a func- 
tion of Vo/Er, for the particular case of *^Rb atoms, A'' — 
2 X 10^, = 27r X 4Hz, lo± =2n x 40Hz, and d = O.S/^m. 



properties and the excitation spectrum are significantly 
different than the expected results for a harmonically 
confined BEC. 
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tial cannot be approximated by the LLL chemical poten- 
tial. As a consequence of that, the challenging problem 
of describing the J//i ~ 1 regime demands a completely 
different approach. In this sense, especially interesting 
could be to employ the analogies to the problem of cou- 
pled Luttinger liquids 2(1] . 

The existence of the quasi Tonks regime can be eas- 
ily revealed in experiments, by either observing the size 
and the form of the stationary density profile, or moni- 
toring the collective excitations. For increasing values of 
Vo/Eji, the frequency of the breathing mode approaches 
that of the lowest compressional mode in ID gases 8], 
whereas the quadrupolar frequencies tend to zero. For 
intermediate values of the lattice potential, the lowest 
excitations are significantly different than in ID, for the 
same value of the dimensionless density 77. In particular, 
more than one excitation mode is expected. We stress 
that the modes (|17ll9f) also differ significantly from the 
expected results for a BEC in a lattice ,18|. The latter 
are recovered from our formalism for large values of A 
(MF regime). We must point, however, that the local 
correlation properties are preserved in the quasi Tonks 
regime. In this sense, if the parameter A becomes suf- 
ficiently small the lifetime of the gas can be very signif- 
icantly enlarged, due to the reduction of the two- and 
three-body losses [lo). 

In this Letter, we have considered what we call the 
quasi- Tonks regime, in which a gas confined in a 2D op- 
tical lattice, can present significant tunneling, and at the 
same time maintain the local chemical potential obtained 
for each site using the LLL model. We have discussed 
that this regime should be the one encountered under 
current typical experimental conditions. In this cross- 
dimensional regime, a 3D cloud acquires ID properties. 
In particular, we have shown that both the ground-state 
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